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Abstract. - We investigate the effect of a rotating Bose-Einstein condensate on a system of 
immersed impurity atoms trapped by an optical lattice. We analytically show that for a one- 
dimensional, ring-shaped setup the coupling of the impurities to the Bogoliubov phonons of the 
condensate leads to a non-trivial phase in the impurity hopping. The presence of this phase can 
be tested by observing a drift in the transport properties of the impurities. These results are 
quantitatively confirmed by a numerically exact simulation of a two-mode Bose-Hubbard model. 
We also give analytical expressions for the occurring phase terms for a two-dimensional setup. 
The phase realises an artificial magnetic field and can for instance be used for the simulation of 
the quantum Hall effect using atoms in an optical lattice. 
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There are several phenomena in physics such as the 
quantum Hall effect [1] or high-temperature supercon- 
ductivity [2], that have been known for more than two 
decades, but where a full theoretical understanding of the 
underlying mechanisms is still lacking. This is due to the 
fact that a complicated many-body system consisting of 
electrons in a crystal has to be described on a quantum 
mechanical level, which quickly exceeds the capability of 
modern computers. As Feynman suggested [3] it is there- 
fore worthwhile to search for alternative quantum systems 
that allow for a clean realisation of the underlying models 
and that are easy to manipulate. With such a quantum 
simulator, we would be able to do simulations that by far 
exceed what is currently possible on classical computers. 

Ultracold atoms in optical lattices are a top candidate 
when it comes to the simulation of condensed matter phe- 
nomena [4]. It has been shown theoretically that for a 
suitable experimental setup the atoms in the lattice in- 
deed realise typical quantum Hall states [5-7]. In a solid- 
state setup, these states are realised by the application of 
a strong magnetic field. Since the motional state of the 
neutral atoms in the optical lattice does not couple to a 
magnetic field, the effect of this field has to be simulated 
by other means. Proposals include using special lattice 
setups and Raman assisted hopping [8], exploiting light 
with orbital angular momentum [9,10], or simply rotat- 
ing the lattice [11]. Whereas the first proposals are rather 
complicated to implement, the rotation of the lattice has 



already been demonstrated in experiments [12]. However, 
in order to observe the quantum Hall states the centrifugal 
term caused by the rotation has to be carefully balanced 
by an additional harmonic trapping potential confining the 
gas, which is difficult to achieve in an experiment. 

In this work, we therefore propose a way of exploiting 
a rotating Bose-Einstein condensate (BEC) to create an 
artificial magnetic field in an optical lattice system which 
is submerged into the condensate. This has the advan- 
tage that the system in which the quantum Hall states 
are to be observed, i.e. the optical lattice, does not need 
to be rotated and thus there is no need of employing a 
balancing potential. The BEC, on the other hand, only 
provides the artificial magnetic field and we are not inter- 
ested in a direct observation of its states. Hence, we can 
employ a harmonic trapping potential which overcompen- 
sates the centrifugal force instead of exactly canceling it, 
which would not be possible if we wanted to observe the 
quantum Hall states in the BEC directly. 

The presence of the condensate can not only be used to 
provide an artificial magnetic field, it also allows for the 
simulation of other effects which are present in condensed 
matter physics [13]. For example, it has been shown that 
for suitably chosen parameters the coupling to the BEC 
phonons may change the transport behaviour of the lattice 
atoms from coherent to incoherent [14] and that an attrac- 
tive interaction between the lattice atoms can be mediated 
by the condensate which leads to clustering effects [15]. By 
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including an artificial magnetic field, the whole setup can 
therefore be used as a well-suited simulator for phenomena 
encountered in condensed matter physics. 

Theoretical description. — In the following, we will 
derive an effective Hamiltonian for the atoms in the opti- 
cal lattice. A similar derivation has been given in ref. [13], 
where it was assumed that the wave function describing 
the condensate is real. In general, this will no longer be 
the case if the condensate is rotated or translated, lead- 
ing to an additional phase in the hopping terms of the 
lattice atoms. The Hamiltonian describing the whole sys- 
tem consists of three parts, H = Hb + Hj + H a . Here, 
H a determines the free dynamics of the atoms of species 
a trapped in the optical lattice, which we will call impu- 
rities in the following. In the laboratory frame, the BEC 
Hamiltonian Hb describing the atoms of species b and the 
density-density interaction Hamiltonian Hj are 

H B = fdr $ ( r ) [#„ + l$( (r) (r) ] ^ (r) ; (1) 



(2) 



where x(r) is the impurity field operator, 4>(r) is the con- 
densate atom field operator satisfying the usual bosonic 
commutation relations, and Ho contains the kinetic energy 
term — fi, 2 V 2 /2m;,, an external trapping potential 14 x t(r), 
the chemical potential fi and any other terms that cor- 
respond to rotations or translations of the BEC as will 
be detailed later. The coupling constants g > and k ac- 
count for the boson-boson and impurity-boson interaction, 
respectively, and m is the mass of a condensate atom. 

In the tight-binding approximation, the impurity field 
operator can be expanded as x(r) = Ej'Tj'M^j) where 
ctj creates an impurity atom in lattice site j and r/j is 
the corresponding Wannier function [16]. The probabil- 
ity densities of these functions have a negligible overlap, 
i.e., J |?7j(r)| 2 |77j/(r)| 2 dr ss for j ^ j 1 , which will be im- 
portant later. By using the above expansion for the field 
operator x the Hamiltonian for the impurities is given by 
H a = -J a J2(i,j) ajaj + (U a /2) £V hj(hj - 1) - fi a E 3 
where \i describes the chemical potential, J a is the hop- 
ping matrix element between two neighbouring sites of the 
lattice, and U a is the on-site interaction strength. 

To investigate the influence of the rotating BEC on the 
impurities we employ the Bogoliubov approximation to 
simplify the dynamics of the condensate. To this end, we 
define </>(r) = </>o(r) + S<fi(r), where <po is a solution of the 
Gross-Pitaevskii equation (GPE) [17] for k = 0, 



= (H + g\M*)\ 2 )Mr) ■ 



(3) 



dimensions. In this case, the deviation of 4>(r) from 4>o( r ) 
is of order k, i.e., (S(j)(r)) cx k, where ( ■) stands for the 
expectation value. Inserting </>o(r) + 5(f>(r) into the Hamil- 
tonian Hb + Hi and keeping terms up to second order in 
k, we obtain the linear term n Jdr x ( r )x( r ) [<^o (r)5<pT (r) + 
05 (r)(50(r)] , in addition to the standard constant and 
quadratic terms in <5</>(r) and S(fv(r). 

The quadratic terms in S(f>(r) are diagonalised by mak- 
ing use of the Bogoliubov transformation. The perturba- 
tion 6<fr(r) is expanded in terms of Bogoliubov phonons 
u q and v q , 5<j>(r) = J2' q l u g( r )b q ~ v*{r)b q ], where b\ cre- 
ates a Bogoliubov phonon in mode q, the prime at the 
sum indicates that the mode corresponding to the ground 
state is not included, and the functions u q and v q solve 
the Bogoliubov-de Gennes (BdG) equations 



Ho + . 9 |^(r)| 2 j «,(r) - g (0(r)) 2 v q (v) = E q u q {v) , (4) 
Hi + g\H*)\ 2 ] Mr) - 9 (0* (r)) 2 u q {v) = -E q v q (r) . 



(5) 



Putting this expansion for the condensate field operator 
into the interaction Hamiltonian and using the fact that 
the overlaps of two different impurity modes is negligible 
allows us to rewrite the total Hamiltonian as a Hubbard- 
Holstein model [18, 19] 



j i 



(6) 



with hw q = E q the energies of the Bogoliubov 

elements Mj q = 
and the 



excitations, the number operator 
the dimensionless matrix 
(n/hw q )Jdr [(f>*(r)u q (r) - </> (r)v 9 (r)] |^(r)| 
mean field shift Ej = k J dr no(r)\r)j (r)| 2 . 

The total Hamiltonian can be brought into a more 
intuitive form by applying the unitary Lang-Firsov 
transformation [18] H c s = UHU\ with U = 
exp [Ei E' q ( M iA ~ Mjjjfy] , which yields 



H eS =UH a fr + J2(Ej - Ejfc - £ Ejh^nj - 1) 

J 3 



(I) 



In contrast to [13] we here do not assume that the wave 
function <f>o is real. To proceed we require a weak impurity- 
boson coupling k such that \K\/gno(r)^(r) <C 1, with 
£h(r) = h/ ^m b gno(v) the healing length, n (r) = |0 o (r)| 2 
the density of the condensate and D the number of spatial The transformation of the impurity Hamiltonian H a 



The presence of the BEC mediates a non-retarded interac- 
tion V h y = Eg huj q (M jtq Mp >q + M^Mj, >q ) between im- 
purities in different lattice sites j and j'. The charac- 
teristic potential energy of an impurity in the deformed 
condensate is described by the polaronic level shift Ej = 
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yields 



UH a V^ = -J a ^2(X i a i ^X j 



U a 



(8) 



where X.] = exp [J2 q (. M j, q b\ ~ M j,qb q )] is a Glauber dis- 
placement operator that creates a coherent phonon cloud 
around the impurity. Thus, the Hamiltonian H e s de- 
scribes the behaviour of polarons according to an ex- 
tended Hubbard model [4] provided that c 3> aJ a /h, with 
c ~ y/gno/rrib the phonon velocity and a the lattice spac- 
ing [13]. 

In ref. [13] it has been shown that for low tempera- 
tures ksT <C Ej the behaviour of the polarons is essen- 
tially coherent. Here, fee is Boltzmann's constant. If also 
J a <C Ej, we can treat the hopping term in eq. (8) as 
a perturbation and trace over the condensate degrees of 
freedom [13]. This allows us to derive the Hamiltonian 



= ~Ji,j a\«j + Uibjfrj - 1) 



(9) 



if- j 



with Jij = fi — kuq + Ej, Uj = U a — 2Ej, and Jjj = 
J a ((xjXj)). Here (( ■ )) denotes the average over the ther- 
mal phonon distribution. After some tedious algebra we 
find 

((XjXj)) = exp \Mi, g - Mj, q \ 2 (2N q {T) + l)j 

x ex P (l Y! M ^ M M ~ M ti M ^J > ( 10 ) 

where N q (T) = [exp(fruj q /kBT) — 1] _1 is the number oc- 
cupation of mode q. As was discussed in ref. [13], the first 
exponential term leads to a suppression of the coherent 
hopping J a . The second exponential, which does not occur 
for the cases presented in [13] since there </>o was assumed 
to be real, has always an absolute value of 1 and causes a 
phase whenever an atom hops from one lattice site j to a 
neighbouring one i. By defining the reduced hopping con- 



stant J a = J a exp 



Mi 



\2N q (T) + l)/2 



the hopping term of the Hamiltonian H^ 1 ' can be writ- 
ten in the form — J a Y2u j) exp(27riaj )J ')a|aj with aiij = 
l/(47ri) J2' q (M itq Ml q - M* q M jtq ). 

This phase depends on the properties of the order pa- 
rameter 4>q and drastically affects the dynamics of the im- 
purities. For example, for a quasi one-dimensional setup 
the presence of such a phase twist in nearest neighbour 
hopping corresponds to a rotation of the ring and thus 
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Fig. 1: (a) Setup of the scheme. The ring confining the BEC 
is rotated with angular velocity Q. This induces a phase twist 
a a when an impurity hops from one lattice site (blue dots) to a 
neighbouring one. (b) Induced phase twist a a versus the angu- 
lar rotation speed Q and the interaction within the BEC. We 
assumed a BEC of 87 Rb with linear density no = 5 x 10 6 m _1 
in a ring with a circumference of L — 12um. The impurities 
are Na atoms in a lattice with 30 sites, which corresponds to 
a lattice constant a = 400nm. The coupling strength between 
the impurity and the BEC is given by k,/2o,Er = 0.035, where 
Eji = (2nh) 2 /2m a (2a) 2 is the recoil energy with m a the mass 
of the impurities. 



causes persistent currents [20]. This is similar to the be- 
haviour of an electron in a superconducting ring subject 
to a magnetic flux. In a two-dimensional setup a suitable 
phase factor realises an artificial magnetic field, which en- 
ables for instance the investigation of quantum Hall ef- 
fects [6]. In the remainder of this paper, we will charac- 
terise the occurring phase term for different setups and 
discuss ways of detecting its influence on the dynamics of 
the impurity atoms. 

BEC in a ring. — We first consider the case where 
the BEC is loaded into a quasi one-dimensional trap of 
length L with periodic boundary conditions, which cor- 
responds to a ring of radius R = L/2ir and has already 
been achieved in experiments [21,22]. The impurities 
are correspondingly trapped by a one-dimensional opti- 
cal lattice which is ring-shaped, see fig. 1(a). We assume 
that the ring trapping the condensate is rotated with an 
angular speed f2. The Hamiltonian Hq describing the 
interaction-free part of the BEC is then given by Ho = 
-(h 2 /2m b )d 2 /dx 2 +ihvd/dx-fx, where v = RCl. The full 
GPE is solved by the wave function 4> = ^/no cxp(iqox) , 
where go — 2irj/L for some integer j such that go = 
rribv/h — Ag with Ag 6 [—tt/L,tt/L]. This choice of go 
ensures that the BEC is in its ground state with a chemi- 
cal potential ^ = h 2 q 2 /2m b +gno—hvqo and that the phase 
of the BEC fulfills the periodic boundary conditions. 

The BdG equations are solved by the mode functions 
u q {x) = A 9 exp[i(g + q Q )x]/y/L and v q (x) = B 9 exp[i(g - 
qo)x]/y/L, where the coefficients are given by the relation 
A q ±B q = (E q /e°) ±:L / 2 , and the quasi-momenta fulfil the 
condition g = 2nj/L for some integer j. The energies of 
the Bogoliubov modes are given by frio q = E q —h 2 qAq/mb 



with E^ = ^e° q (e° q + 2gn ) and e° = h 2 q 2 /2m b . In or- 
der to calculate the coupling matrix elements Mj. q we as- 
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sume that the lattice trapping the impurities is sufficiently 
deep such that the Wannier functions rjj are well approx- 
imated by Gaussians of width a centered at lattice site j, 
i.e., r)j(x) ~ exp[— (x — Xj) 2 /2a 2 ]/ \J y/ira, where Xj is the 
position of lattice site j. We furthermore note that due 
to the translational invariance the phase a%i+i — a a does 
not depend on the index i and we thus find 

a a = - — > % e q i sm(ga) . 11 

q^qo v 

We see that the finite width of the Wannier functions a 
leads to an effective cut-off of the sum for large quasi- 
momcnta q and thus large energies hw q . Furthermore, 
the expression for a a contains the Bogolibov excitation 
energies hio q of the rotating condensate as well as the static 
structure factor e q /E q of the non-rotating BEC [23,24]. 

The value of a a depends critically on the rotation speed 
Q. For zero rotation it can be shown from eq. (11) that 
a a = 0. If the rotation speed is increased, the phase 
twist a a initially increases linearly with S7, as shown in 
fig. 1(b). However, for a critical rotation speed a jump 
occurs and the value of the phase twist changes from pos- 
itive to negative. This behaviour can be explained as 
follows: Due to the finite size of the ring, the possible 
quasi-momentum states of the BEC are restricted to values 
qo = 2nj/L for some integer j. For rotation speeds close 
to zero, the ground state of the BEC exhibits zero quasi- 
momentum. Only if the rotation speed is above a critical 
value f2 C rit = h/2m,bR 2 will the ground state change to 
a non-zero quasi-momentum. Further quasi-momentum 
jumps occur at odd multiples of this critical angular speed. 
Together with the change in the quasi-momentum of the 
BEC, also the jump in the induced phase twist a a occurs. 
To some extend one could thus say that the induced phase 
is caused by a mismatch of the angular rotation speed 
and the quasi-momentum q of the BEC, which is enforced 
to be quantised by the boundary conditions. It should be 
noted that at the critical points where the BEC changes 
its quasi-momentum, the ground state of the condensate 
is degenerate, which leads to a failure of the Bogoliubov 
approximation [25]. 

The parameters presented in fig. 1(b) fulfill the condi- 
tions for our derivation to be valid and can be achieved 
using present experimental techniques [13-15,26,27]. The 
precise value of the coupling constant g can be changed 
via a Feshbach resonance. The temperature chosen for 
our calculations is zero, but the discussed effects will be 
observable as long as the temperature is lower than the po- 
laron energy, which is typically on the order of a few tens 
of nanokclvin [13]. This also ensures that the first factor 
in eq. (10) is close to one leading to sufficiently large hop- 
ping. In the above examples the size of the induced phase 
twist can reach values as large as a a ~ 0.03. This will 
significantly affect the phase of the atomic wave function 
and make effects of the induced artificial magnetic field 
observable even if the lattice only consists of a few tens of 



sites in each direction. 

Drift of the impurities. The presence of the phase 
twist a a can be measured by observing its influence on 
the transport properties of the impurities. We consider the 
following case: A single impurity is initially trapped in a 
finite spatial region of the lattice. This can be achieved by 
applying an additional trapping potential Vi n i(j) = Vo(j — 
jo) 2 , in which case the single impurity starts in the state 
l"0)ini oc £j ex P [-{j - j ) 2 /2w 2 ni ] a]|vac), where w mi is 
the width of the initial state. Without a phase twist a a 
after turning the initial potential off, the impurity expands 
symmetrically around the lattice site jo and no overall 
drift to either higher or lower lattice sites occurs. 

This behaviour changes if a finite phase twist a a is taken 
into account as shown in fig. 2(a). For the chosen pa- 
rameters and a relatively short time of t — 5h/ J a , which 
under standard experimental conditions corresponds to a 
few miliseconds, a drift towards higher lattice site num- 
bers is clearly visible. Such drifts can be measured in 
experiments [21] and thus allow to probe the presence of 
an induced phase twist. The mean position of the im- 
purity after an evolution time of t = 5h/J a is shown in 
fig. 2(b). It exhibits a periodic behaviour in a a with a 
period of 1, which is caused simply by the fact that the 
integer part of a a only contributes a trivial phase. In to- 
tal, the behaviour is akin to an electron in a ring subject 
to a magnetic flux. The electron will also drift towards 
a preferential direction and thus will cause a current as 
known from superconducting rings. 

Numerical simulations of the full ring system. We now 
compare our prediction to a system where a numerically 
exact solution is possible, namely a small two-mode Bose- 
Hubbard model. We assume that the condensate is rep- 
resented by a mode c that experiences a rotation whereas 
the other mode a representing a single impurity docs not. 
The system is described by the Hamiltonian 

H B h =~ a\a,j +Ui^2 aj&jcjcj- (12) 

-J C J2 (e 27riQ ^+i + h.c) + ^ Yl ¥i c & ■ 

J 3 

Here, J c is the hopping matrix clement of the atoms in 
mode c, U c their interaction constant, Ui the interaction 
between atoms in mode c and a, and a c = m c n c N s a 2 /2irh, 
where m c is the mass of the atoms, f2 c the speed with 
which the ring for atoms in mode c is rotated, N s the 
number of lattice sites, a the lattice constant, and periodic 
boundary conditions have been assumed [25]. 

We first calculate the ground state of the atoms in mode 
c for different values of a c and vanishing coupling to the 
impurity atom. We assume that the impurity is in state 
\ip) ini at time t = when the interaction {//is turned on 
and let the system evolve for a finite evolution time. For 
two different fl c cx a c the density distributions of the im- 
purity after this time arc shown in fig. 3(a), using typical 
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Fig. 2: (a) Plot of the impurity distribution versus dimension- 
less time td = t.Ja/h. A darker colour corresponds to a higher 
density, (b) Mean position of the impurity after an evolution 
time of td = 6. We have chosen an initial state \tp)i a i of width 
Wini = v2 centered at lattice site jo = 15. In (a), the chosen 
phase twist was a a = 0.03. 

experimental parameters. A small drift to either the left 
or the right is clearly visible. We expect this difference 
to be by far more pronounced for a larger lattice with a 
longer evolution time. However, for computational rea- 
sons we had to restrict ourselves to a small system size. 
The behaviour of the mean end position versus the phase 
twist a c is shown in fig. 3(b). As we see this drift exhibits 
a jump at a value of a c = l/2N s . The reason for this 
jump is similar to the saw-tooth behaviour of the induced 
phase a a observed in fig. 1 (b) . Due to the quantisation of 
the BEC quasi-momenta there exist critical rotation ve- 
locities fi£ rit = 2n 2 (2j + l)h/m b N 2 a 2 at which the ground 
state of the BEC changes from one quasi-momentum state 
to another one [25]. These critical velocities correspond 
to phase twists a c = (2j + 1)/2N S for integer j. One 
can show that due to these jumps the induced phase a a 
changes analogously to the behaviour shown in fig. 1(b). 
This change is then also reflected in the drifts that the 
impurity atoms experience. 

The two-dimensional condensate. A phase twist 
a»,j is also induced if we consider a two-dimensional BEC 
in a trap rotating with angular velocity Q. For simplic- 
ity we assume an interaction free condensate, i.e., g = 0. 
In the rotating frame, the (single particle) Hamiltonian 
describing the condensate subject to a harmonic trap- 
ping potential u> tT thus reads Hq = — (ft 2 /2m;,)(V') 2 + 
mf,w 2 r (r') 2 /2 — fj, — fle z • r' x p', where e 2 is the unit 
vector in z direction and p' = — ihV. It can be shown 
that the solutions for this Hamiltonian arc given by 
i> n ,l(r,<p) = y/n\/n(n + 1)1 e^e"^ ^ L 1 ^ 2 ) I '<*) with en- 
ergies E n> i+fj, = {2n+l+l)fajJt, t —ttll. Here, £ = r/a with 
ao = ^Jh/mfjUJtx the harmonic oscillator length and (r, ip) 
are polar coordinates. The quantum numbers n and I are 
restricted to n = 0, 1, 2, ... and I = — n, —n + 1, ... . From 
this we see that for high rotation speeds |fi| > u> tT the 
system is unstable since for suitable (n, I) the energy can 
be made infinitely small which corresponds to the atoms 
being driven infinitely far away from the origin due to the 
centrifugal term. We thus restrict the rotation speed to 
moderate values |f2| < w tr , for which the ground state is 




Fig. 3: Expansion of a single impurity atom initially in state 
IV^mi with jo = 11 and width uiini = V2. (a) Density distri- 
bution after a time td = tj a /h — 4 for q c = 0.0227 (dashed 
line) and a c = 0.0251 (solid line), (b) Mean position of the im- 
purity after the evolution time td = 4. The two crosses mark 
the positions of the densities in (a) . The other parameters are 
U c = Jc = Ja, Ui = 2J C with three atoms in mode c and a 
lattice consisting of N 3 = 21 sites. 

always given by (n, I) = (0, 0). 

For brevity, let us introduce the multi- index v = (n,l). 
If we assume that the condensate is in a state vq, we can 
expand the field operator of the condensate as ip(r, ip) — > 
VNo^u {r,ip) + E^o^vfov)^, where (Sj, O = N 
is the number of atoms in state vq, and we assume that 
{b\,b u ) <C A^o for v ^ vq. As we consider the interaction- 
free case this is exactly the same as an expansion in terms 
of Bogolibov phonons, which would yield v v = and 
iij, = ip v (r, ip) . The chemical potential is determined by 
solving the Gross-Pitacvskii equation for the correspond- 
ing state vq, which gives [i Va = —E n[u i Q . This allows us to 
derive the phase twist ajj in the same way as for the ring 
optical lattice. In order to solve the integral in the ex- 
pression for the coupling matrix elements Mj^ we assume 
that the Wannier functions are strongly localised and can 
be approximated by a Dirac delta function. If the con- 
densate is in the ground state v$ = (0,0), we find for the 
induced phase twist 

(o,o) _ k 2 N ^ 1 n\ -t?-^ 

-2^at^ o (h^) 2 (n + l)! C (13) 

where (Ci > Vi) and (£j,(/2j) are the coordinates of the two 
lattice sites before and after the jump of the impurity. We 
immediately see that if ifi = (fj, which corresponds to 
a hopping in radial direction, the induced phase is zero. 
Furthermore, the induced phase vanishes as £j or £j go to 
infinity, which can be explained by the vanishing density 
of the condensate for large distances. 

Our theory is also applicable if the condensate is initially 
in an excited state, for example in the state vq = (0, 1) 
that exhibits a vortex. In this case, we derive for the 
phase twist 

(o,i) _jfNo_ v _L 2!_ e -c?-cj^. f .y+i 

' 3 ^HX(^) 2 (n + iy (14) 
xLL(4 2 )^(^)sin[(^-^)a + l)]- 
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Fig. 4: (a) Reduced phase twists = a ij x 

(K 2 N /2Tv 3 at(hujtv) 2 )~ 1 for the ground state (n,l) = (0,0) 
(solid line) and the vortex state (n, I) = (0, 1) (dashed line) 
versus dimensionless angular rotation speed fl/c^tr- Here, 
£i = £j = 0.5 and <fii — (pj — 0.1. (b) Reduced phase twists ver- 
sus £ = = fj- The angular velocity is given by fi = 0.25tJtr 
and tpi — ipj — 0.1. 

The behaviour of the two phase twists for varying an- 
gular rotation speed is shown in fig. 4(a). In both cases, 
the induced phase twists increase as the rotation speed 
£1 is increased and diverge as £1 approaches the trapping 
frequency w t r ■ We furthermore observe that for small rota- 
tion frequencies the state exhibiting a vortex still induces 
a finite phase twist, whereas oif- ^ vanishes as f2 goes to 
zero. The dependence of the phase twist on the distance 
from the origin is shown in fig. 4(b). For large distances 
£, the phase twist is zero due to the vanishing conden- 
sate density. For small distances the phase twist vanishes 
as well, due to the fixed difference ifi — ipj. The spatial 
dependence of a£°- leads to the realisation of an artificial 
magnetic field, which in general will be non-homogenous. 
We expect that the spatial dependence of this field can be 
controlled by varying the trapping potential and thus the 
density of the BEC. The behaviour of the impurities can 
also be used to probe the state of the condensate. These 
effects will be the topic of future research. 

Summary. — We have investigated an alternative way 
to create artificial magnetic fields in a system of ultracold 
atoms in an optical lattice. We have shown analytically 
that for impurities trapped in a ring-shaped optical lattice 
submerged into a ring-shaped, rotating Bose-Einstein con- 
densate a drift in the transport of the impurities occurs. 
This drift was qualitatively confirmed by a full numerical 
simulation of a two-mode Bosc-Hubbard model. We also 
showed that an artificial magnetic field can be created in 
a two-dimensional optical lattice by submerging it into 
a moderately rotating, harmonically trapped condensate. 
With this, our system lends itself as a well-suited quan- 
tum simulator for investigating phenomena encountered 
in condensed matter physics. 
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